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Abstract. Contrary to the common lore based on naive dimensional analysis, the seesaw 
scale for neutrino masses can be naturally in the TeV range, with small parameters coming from 
radiative corrections. We present one such class of type-I seesaw models, based on the left-right 
gauge group SU{2)l x SU{2)r x U{1)b-l realized at the TeV scale, which fits the observed 
neutrino oscillation parameters as well as other low energy constraints. We discuss how the 
small parameters of this scenario can arise naturally from one loop effects. The neutrino fits 
in this model use quasi-degenerate heavy Majorana neutrinos, as also required to explain the 
matter-antimatter asymmetry in our Universe via resonant leptogenesis mechanism. We discuss 
the constraints implied by the dynamics of this mechanism on the mass of the right-handed 
gauge boson in this class of models with enhanced neutrino Yukawa couplings compared to the 
canonical seesaw model and find a lower bound of mwfi > 9.9 TeV for successful leptogenesis 
assuming maximal CP asymmetry for each flavor. We also present a model with explicit neutrino 
mass fit, where the lower bound goes up to 13.1 TeV due to less than maximal primordial CP 
asymmetry predicted by the model. 


1. Introduction 

The seesaw mechanism [T] seems to provide a simple way to understand the smallness of neutrino 
masses. The simplest among the various seesaw mechanisms is the type-I seesaw, which requires 
the introduction of Standard Model (SM) gauge-singlet right-handed (RH) neutrinos (iV’s). The 
Majorana mass Mjy of these RH neutrinos explicitly breaks the global B — L symmetry of the 
SM by two units and gives rise to a light neutrino mass matrix of the form 

M, ^ -MdM^^mI , (1) 

where is the Dirac mass generated due to the neutrino Yukawa couplings after the 
electroweak symmetry breaking. If one assumes maximal allowed values for the third-generation 
Dirac mass as is implied by certain Grand Unified Theories such as SO(IO), then the 

seesaw scale is pushed up to ~ 10^^ GeV making it virtually inaccessible to both collider 
experiments at the energy frontier and other low energy searches at the intensity frontier. 
Such naive (“dimensional counting”) considerations however can be misleading in the context 
of specific ultraviolet (UV) complete seesaw models where there may be dynamical suppressions 
of neutrino Dirac masses. One class of such models that are based on the SM gauge group 
generally go by the name of radiative seesaw models; see e.g. [2]. Within the type-I seesaw 



framework, there exists another possibility to suppress the light neutrino masses by assigning 
specihc textures to the Dirac and Majorana mass matrices in the seesaw formula O 0]. The 
stability of these textures can in principle be guaranteed by enforcing some symmetries in the 
lepton sector. 

Extending the SM gauge group to the Left-Right (L-R) Symmetric gauge group SU{2)l x 
SU{2)ii X U{1)b-l provides a simple UV-complete seesaw model [5]. Here, the RH neutrinos 
are a necessary part of the model and do not have to be added ‘by hand’ just to implement 
the seesaw mechanism. An important point is that the RH neutrinos acquire a Majorana mass 
as soon as the SU{2)fi symmetry is broken at a scale vr. This is quite analogous to the way 
the charged fermions get masses in the SM by the Higgs mechanism when the SU{2)l gauge 
symmetry is broken at a scale v. The Higgs field that gives mass to the RH neutrinos becomes 
the analogue of the 125 GeV Higgs boson discovered at the LHC. 

In this proceedings, we focus on a new class of L-R seesaw models [6], where in addition to 
the RH gauge bosons as well as other Higgs fields, light (TeV scale) RH neutrinos with sizable 
mixing with the left-handed (LH) neutrinos may be accessible at colliders and other low energy 
experiments (for a review, see e.g. [7]). With the LHC ramping up to a higher center of mass 
energy ^/s = 13 TeV and later to 14 TeV, our model as well as other low scale seesaw models, 
should provide motivation for testing the origin of neutrino mass in the coming years. For the 
collider phenomenology of TeV-scale LRSM, see e.g. [EIEIIIQ], and for the current experimental 
status, see m- Such considerations should also be of interest m for setting physics goals for 
future higher energy colliders such as the ^/s = 80 or 100 TeV VLHC |13j . It is also important 
to mention here that the new class of L-R seesaw models with relatively large left-right neutrino 
mixing can be distinguished from the minimal LRSM using different kinematic variables [iniiiij. 

An attractive consequence of the seesaw mechanism is the possibility that it could solve 
another important puzzle of cosmology, i.e. the origin of matter-antimatter asymmetry in our 
Universe, via a mechanism called leptogenesis m- At the heart of this mechanism is the out- 
of-equilibrium decay of the RH Majorana neutrinos via the decay modes N —)• (where 

the superscript c denotes the CP-conjugate), which violate L, C and CP, thereby dynamically 
generating an asymmetry in the lepton sector. This primordial lepton asymmetry undergoes 
thermodynamic evolution as the Universe expands and finally gets converted to a baryon 
asymmetry via equilibrated electroweak sphaleron interactions m- The attractive aspect of 
this mechanism is that the same Yukawa couplings that give rise to neutrino masses via 
the seesaw mechanism in Eq. 0 are also responsible for the origin of matter, thus implying 
an intimate connection between two seemingly disparate pieces of evidence for beyond the SM 
physics. This beautiful idea can be tested, provided the seesaw scale, and hence, the scale at 
which leptogenesis takes place, is accessible to current and near future laboratory experiments. 
Low-scale leptogenesis is possible when at least two RH Majorana neutrinos have a small 
mass difference comparable to their decay widths, thus resonantly enhancing the e-type CP 
asymmetry [HIM due to RH neutrino self-energy effects m- 

An important question to explore is whether it is possible to falsify the leptogenesis 
mechanism by laboratory experiments |20j . In particular, it is worth investigating if the 
leptogenesis mechanism works in a TeV-scale realization of L-R seesaw and whether an 
observation of a RH gauge boson at the LHC could falsify this scenario. This problem was 
analyzed in detail in [21] and it was pointed out that the additional dilution and washout effects 
on the lepton asymmetry due to kUR-mediated AL = 1 scattering processes become important 
for a TeV-scale L-R seesaw. Since in generic versions of the seesaw model, sub-eV neutrino 
masses would require the Dirac Yukawa couplings to be the washout effects would 


lead to an efficiency factor n 


9r^n 


{gji being the SU{2)r gauge coupling, which is assumed 


to be equal to the SU{2)l gauge coupling in the minimal scenario) which is too small for a low- 


2 



scale leptogenesis to work. Through a detailed analysis of the relevant Boltzmann equations, it 
was concluded in [21] that even for maximal CP-asymmetry e ~ 0(1), the observed value of the 
baryon asymmetry can be explained by leptogenesis in these L-R models, only if > 18 TeV. 
Turning this argument around, if a positive signal for Wr is observed at the LHC, this will falsify 
leptogenesis as a mechanism for understanding the origin of matter in L-R seesaw framework. 
Since this is such an important issue, we have reanalyzed this |22j to examine the robustness 
of the lower bound on and to investigate if there exists any allowed parameter space with 
successful leptogenesis for smaller values of • We find that for the generic class of L-R seesaw 
models with large light-heavy neutrino mixing, successful leptogenesis requires mw^i > 9-9 TeV 
for the maximal CP asymmetry. For illustration, we will present an explicit neutrino mass fit 
which satisfies the observed baryon asymmetry. The main new results presented here should 
supersede those reported in [22] . 

This proceedings is organized as follows: in Section we review the basic features of the 
generic L-R seesaw models as well as our new L-R seesaw model with special Dirac and Majorana 
textures resulting in light neutrinos via type-I seesaw with large light-heavy neutrino mixing 
generated in a natural manner. In Section]^ we present an alternative neutrino mass fit relying 
on cancellations in the seesaw formula. In Section [^ we discuss resonant leptogenesis in this 
model and the impact of a TeV-scale Wr on the lepton asymmetry. Our conclusions are given 
in Section m 


2. A TeV-scale left-right seesaw model 

In the minimal L-R symmetric model (LRSM), the fermions are assigned to the gauge group 
SU{2)l X SU{2)r X U{1)b-l as follows: denoting Qi = {u,d)l and V'i = as the quark 

and lepton doublets of the zth generation respectively, and (also denoted simply by 
Li) are assigned to doublets under the SU{2)r group, while Qr^i and ^pR^i (also denoted by Ri) 
as the doublets under the SU{2)r group. The Higgs sector of the model can consist of one or 
several bidoublets (pa and triplets Ar^i,: 


P = 




Ar 


~^r/^ ) 


( 2 ) 


The gauge symmetry SU{2)r x U{1)b-l is broken by the vacuum expectation value (VEV) 
(A^j) = vr to the group U{1)y of the SM. The LH counterpart (A/,) to Ar is not considered 
here, assuming that parity and SU{2)r gauge symmetry scales are decoupled. In this case, the 
Ar fields become heavy when the discrete parity symmetry is broken, and disappear from the 
low energy theory [23|. The VEV of the p field given by {p) = diag(«:, k') breaks the SM gauge 
group to t/(l)em- The fermion masses are obtained from the following Yukawa Lagrangian: 

-Cy = hjfQRipaQR^j + hff‘QL,i^aQR,j + UPaRj 

+ hpj LipaRj fij{RjCiT 2 ARRj lJC iT 2 ARLj) -|- H.c., (3) 


where p = T 2 p*T 2 (t 2 being the second Pauli matrix). After electroweak symmetry breaking, the 
Dirac fermion masses are given by the generic formula Mf = n + k' for up-type fermions, 
and for down-type quarks and charged leptons, it is the same formula with k k'. The Yukawa 
Lagrangian (j^ leads to the Dirac neutrino mass matrix Mr) = + HJ-k' and the Majorana 

mass matrix for the heavy RH neutrinos Mjy = fvR which go into the seesaw formula ([^ for 
calculating the neutrino masses and the heavy-light neutrino mixing. 

Since the Higgs sector relates the neutrino Yukawa couplings with the charged-lepton ones in 
the LRSM , it is interesting to see how a TeV-scale L-R seesaw can be realized in Nature. There 
are essentially three ways to do so: (i) by choosing one set of Yukawa couplings to be < 


3 


for a particular VEV assignment for the bidoublet Higgs fields; (ii) by choosing particular 
symmetry-motivated textures for the Yukawa couplings which ensure that the leading order 
seesaw contribution to neutrino masses vanish; and (iii) by choosing large Yukawa couplings as 
allowed by specific textures and invoking cancellations in the seesaw formula 0. The strategy 
(i) with small Yukawa couplings does not lead to any new interesting phenomenology over the 
minimal LRSM. Therefore, we will focus on strategies (ii) and (iii) in the following sections, and 
in particular, use strategy (iii) for our numerical analysis for the calculation of lepton asymmetry. 
A similar analysis can be performed following strategy (ii) as well. For the phenomenological 
implications of this new class of L-R seesaw models, see e.g. mm- 


2.1. A symmetry-protected L-R seesaw with large light-heavy neutrino mixing 
As mentioned above, the basic strategy for understanding small neutrino masses with a low-scale 
type-I seesaw is to have the appropriate textures for Mjo and Mjq in Eq. Q. There are several 
examples of this type laii! in the context of the minimal seesaw extensions of the SM. Here, we 
discuss the embedding of one such texture (from Ref. [3] ) in the L-R model using an appropriate 
family symmetry [ 6 ]. The symmetry must not only guarantee the special leptonic textures but 
also must be free of light scalar bosons which can result if the effect of the discrete symmetry is 
to automatically lead to a 17(1) symmetry. Moreover in LRSM, the charged-lepton mass matrix 
Ml is related to Mjo which puts additional constraints on phenomenological viability of the 
model. Therefore, we find it remarkable that the model presented below remains a viable TeV- 
scale L-R type-I seesaw model for neutrinos, and as a result, has interesting phenomenological 
implications beyond the minimal LRSM. 

The Dirac and Majorana mass matrices Mjo and M^r considered here have the following form: 


Md 


/ mi 

<^1 



/ 6M 

Ml 

0 \ 

m2 

^2 

62 , 

Mat = 

Ml 

0 

0 

\ m 3 

<^3 

63 / 


\ 0 

0 

M 2 / 


(4) 


with €i,6i <C mi and 6M <C Mj. In the limit of €i,6i,6M —>■ 0, the neutrino masses vanish, 
although the heavy-light mixing parameters given by ^ij = mi/Mj can be quite large. The 
light neutrino masses given by the seesaw formula a become proportional to either e*, 5i or 
their products. If the smallness of 5i and can be guaranteed by some symmetry, then we 
have a “natural” TeV-scale seesaw model with tiny neutrino masses. As we show below, the 
mass textures Q can be successfully embedded into the L-R framework, while satisfying all 
current experimental constraints. Note that the choice of texture for the RH neutrino mass 
matrix in Eq. Q is necessary to keep the model compatible with strong bounds from tree-level 
IjL~ —)> e“e'''e“ decay via TeV-scale Ar exchange [22] . 

In order to obtain the special Dirac and Majorana textures given in Eq. Q, we supplement 
the L-R gauge group SU{2)l x SU{2)ji x U{1)b-l with a global discrete symmetry D = 
Z 4 X Z 4 X Z 4 [ 6 ]. For the Higgs sector, we choose three bi-doublets (i;i>i, 2 , 3 ) with H — L = 0 and 
two RH triplets Ab, 2 ) with B — L = 2. The fermion and Higgs multiplets are assigned 

the D quantum numbers as shown in Table The leptonic part of the Yukawa Lagrangian Q 
invariant under the D-symmetry is given by 

— A,y = hiiLi^iRi hi 2 Li(j) 2 R 2 + hi'iLi<j)^R 2 , + fi 2 RiR 2 AR^i -|- fz?,R‘iR‘i^R ,2 + H.c. (5) 

In the symmetry limit, the VEVs of (^ 2 ,^ have the form {(pa) = diag(0, fi:(j). This is because the 
terms of the form AT{pa(p\) which would change the (p VEV to the form diag(K, n') are forbidden 
by the D-symmetry from appearing in the scalar potential. Thus in the symmetry limit, the 
Mr elements in one column are big and non-zero, and the charged lepton mass matrix has one 
eigenvalue zero which can be identified as the electron flavor. The RH neutrino Majorana mass 
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Field 

Z4 X Z4 X Z4 transformation 

Ti,2,3 

(1, 1, 1) 

Ri 

{-h 1, 1) 

R2 

(1, -b 1) 

R3 

(1, 1, -i) 

4>i 

{-h 1, 1) 

4>2 

(1, b 1) 

4>3 

(1, 1, i) 

^R,l 

(b b 1) 

^R,2 

(1, 1, -1) 


Table 1: The discrete symmetry assignments for the fermion and Higgs fields in our L-R model 
that lead naturally to the special Dirac and Majorana textures given in Eq. 0 M- 


matrix has the form in Eq. 0 with 6M = 0. To make the model realistic, we add small soft 
symmetry-breaking terms of the form 

+ ( 6 ) 

a^b 


to the scalar potential. This will induce the (j)a VEVs of the form {(pa) = diag((5Ka) n'a)^ where 


5Ka oc 


Eb f^lb4 


Ea,b ^'ab'^' 


a'^'b 


+ Ea E 


^R,a 


(7) 


with A, A' respectively being the scalar self-couplings of the bidoublet and triplet fields in the 
scalar potential. Choosing appropriately small we can get very small 6Ka as required to 
satisfy the neutrino oscillation data. 


2.2. Naturalness of small 6k 

It is possible to generate the small parameters 6Ka naturally through loop effects involving the 
Wl — Wr mixing. To explain this, first we note that if in a L-R model, 6k's are set to zero in 
a natural manner, it will lead to zero Wl — Wr mixing. By the same token, when Wl — Wr 
mixing is nonzero, it will induce a nonzero 6k. Thus there is an intimate connection between 
6k and Wl — Wr mixing. In our model discussed above, Jk = 0 is guaranteed by {Z^)^ discrete 
symmetry. However, if bi-doublets are used to generate quark masses, it is not possible to set 
either k or k' to zero, since in that case the symmetry that stabilizes the vanishing of k or k' 
makes it impossible to get realistic quark masses and mixing. Once in the quark sector both 
K and k' are set to nonzero values, they will induce tree level Wl — Wr mixing and hence 
arbitrary 6k ^ 0 at the tree level in the leptonic sector, thus invalidating the naturalness of the 
whole scenario. There is however an alternative way to get quark flavor pattern right without 
introducing bidoublet Higgs fields and thereby avoiding nonzero Wl — Wr mixing at the tree 
level. If we introduce SU{2)l,r singlet vector-like up {U) and down (V) quarks, and only LH 
and RH doublet Higgs fields xl,r with nonzero VEVs which give mass to the quarks via a quark 
seesaw mechanism as in [23], it keeps the 6k’s naturally small. To briefly introduce the quark 
seesaw, we write the Yukawa interactions responsible for fermion masses in this model as 

— C.Y = YuQlXlUr + YciQlXlT^r + {L R) + UlMuUr + VlMrVr + H.c.. (8) 
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Figure 1: One loop graph for Wl — Wr mixing in the quark seesaw model. The (blue) crosses 
in the middle of the two quark propagators are for the vector like quark mass insertions, while 
the other four (black) crosses denote the mixing between the vectorlike quarks and SM quarks. 


where xl,r = r- This leads to generic seesaw type quark mass relations, e.g. 


Mu 




yuVLVR 

2Mi( 


(9) 


The leptonic sector bidoublets do not couple to the quark sector due to the {Z^)^ discrete 
symmetry under which the quark sector fields are singlets. As a result, the Wr — Wr mixing is 
zero also in the quark sector, and hence, 6k is not induced at the tree level. However, Wr — Wr 
mixing is induced at the one loop level via the exchange of top and bottom quarks in the quark 
sector (see Figure 1), which in turn induces leptonic 6k at the two loop level from the gauge 
interaction gRgR<f>^^(p 2 W^W^ in the notation of Eq. (§. This makes 6k finite and small, thus 
keeping the model technically natural. 

An important point to notice is that 6KaS are responsible for the electron mass as well as 
neutrino masses via type-I seesaw. Thus getting a fit to the observed neutrino masses and mixing 
while at the same time keeping electron mass at its desired value is a nontrivial task since in 
the lepton sector, the model has only 14 free parameters (9 Yukawa couplings, 3 Majorana mass 
parameters and two VEVs) to start with. One can always go to a basis where Mr, takes an 
upper-triangular form, in which case the number of free parameters reduces to 11. Out of this 
11 model parameters, we must not only get fits for the three charged lepton masses, the two 
neutrino mass-squared differences and three mixing angles (total of 8 outputs), but must also 
satisfy the unitarity constraints on the new light neutrino mixing matrix as well as constraints 
from rare lepton decays which involve only the RH mixing matrix for charged leptons. It was 
explicitly shown in that this model does indeed provide a fit to all observables while satisfying 
all the necessary constraints. 


3. An alternative nentrino mass fit 

In this section, we discuss neutrino fit in a model with Mr, and Mv texture as given in Eq. @ 
without worrying about its origin from discrete symmetries and use this fit in the subsequent 
section for the discussion of leptogenesis. We consider a minimal scenario with only one bidoublet 
(p and a single right-Higgs triplet to illustrate the fit. Clearly, in this case we will have to 
adjust parameters to get the desired fermion mass textures. 

In this class of L-R models, we can always choose a basis of the LH sector prior to 
SU{2)r X U{1)y breaking such that the Dirac mass matrix Mr, can be written in an upper- 
triangular form without affecting the RH neutrino texture. In this basis, the mass matrices 
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obtained from the Yukawa Lagrangian in Eq. i) are given by 



( 

hnn' + 

huK* 

hi2K' -1- hl2K* 

hisK' -h hi3K* 

II 


\kT-\k 

n'* 

-/.21 

h22K' + h22K* 

/l23«:' -h h23M 


[ 

|k'P-|k 

k'* 

'/l31 


+ h33K* 


( 

hiiK + Jiiik'* 

hi2K -1- hi2K'* 

hisK -h Jiisk'* 

Md = 


0 


/l22'^ + h22K'* 

^23^ + h23K'* 


[ 

0 


0 

hssn + 


( 

6M 

fl2VR, 

1 0 \ 


Mn = 


fl2VR,l 

0 

0 



[ 

0 

0 

2/33'1’K,2 / 



(10) 


( 11 ) 


( 12 ) 


Note that since parity is a broken symmetry in this low energy version, we need not have the 
Yukawa couplings hij to satisfy hermiticity. We choose the parameters of the model, i.e. hij's 
and K, k' such that all the experimental constraints in the lepton sector are satisfied. As an 
example fit, we consider the mass matrices of the following form: 


Ml 


Md 


Mn 


0.00120 -0.0507 -1.41 \ 

0 0.0929 -0.657 GeV, 

0 0 -0.861 / 

0.0676 -1.06 x lO”’^ _;l.06 x lO"'^ \ 

0 1.95 X 10“^ g-i0.557n _4 7g ^ ;^q- 5 QgY^ 

0 0 -6.38 X 10"^ / 

2.59 X 10-® 585 0 \ 

585 0 0 GeV, 

0 0 -585 / 


(13) 

(14) 

(15) 


In addition, we choose the VEVs k = 6.46 GeV and k,' = 173.88 GeV, which satisfy 

|«:^| + with V = 174 GeV. These are all the input parameters for the model. The Mi 

is diagonalized by a bi-unitary transformation: Mi = (V)'^)1M;V)^, where 


V,^ = i 


/ 0.426 

0.436 

0.793 \ 

/ 1.00 

0.00495 

0.000535 \ 

0.232 

-0.900 

0.369 

, Vi^ = i\ 0.00495 

-1.00 

-0.00332 

V -0.875 

-0.0271 

0.484 / 

V 0.000519 

0.00332 

-1.00 / 


(16) 


The light neutrino mass matrix M^, is diagonalized by a unitary transformation: M^, = VjM^V^^ 
where 

/ 0.338i -0.0149 0.941 \ 

W = 0.742i 0.619 -0.257 . (17) 

\ -0.579i 0.785 0.221 ) 


The Pontecorvo-Maki-Nakagawa-Sakata (PMNS) mixing matrix is then given by 


/ 0.823 0.549 0.148 \ / 1 0 0 \ 

VpMNS = (V)^)^K = -0.505 0.585 0.635 0 z 0 , 

\ 0.262 -0.597 0.758 /VO 0 -i ) 


(18) 


The resulting neutrino masses and mixing angles as well as the charged-lepton masses are given 
in Tabled 
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Parameter 

Value 

nie 

0.511 MeV 

rrifj, 

106 MeV 

IJlr 

1.78 GeV 


6.49 X 10-5 eV 

171^2 

1.08 X 10-2 eV 


5.02 X 10-2 eV 

0l2 

33.7° 

6*23 

40.0° 

6*13 

8.51° 

rriNi 

585 GeV 

mN2 

585 GeV 

rriNs 

585 GeV 


Table 2: Masses and mixing angles calculated from the mass matrices given by Eqs. 
The mass splitting between the RH neutrinos is of order of 6M r\j 10“® GeV. 


(13)-(15). 


Since Mi = Hk' + hn* and = hK + hn'*, the Yukawa couplings h, h are easily calculated 
using VEV’s chosen as above, and they are found to be 


h 


h 


(0.15-k 1.98f) X 10-5 

(1.72 -k 2.36f) X 10-^ (4.77 -k 6.57f) x 

10-3 \ 


0 

-(3.14 -k 4.33f) X 10-"* (2.22 -k 3.06f) x 

10-3 , 

(19) 

0 

0 

(2.91 -k4.01f) X 

10-3 / 


(2.69-k2.81f) X 10-^ 

1.08 X 10-5 

3.02 X 10-^ \ 



0 

-1.99 X 10-5 

1.41 X 10-"* . 


(20) 

0 

0 

1.84 X 10-"* / 




The leptonic CP asymmetry generated by the RH neutrino decays is governed by these Yukawa 
couplings. If we assume that the two SM doublets in the LR bi-doublet i.e. 4>i 4>2 are 

approximate mass eigenstates with > mjy > the dominant contribution to the CP 

asymmetry comes solely from the self-energy correction by the (/) 2 -loop to the decay process 
N —>• 4>2Li. Therefore, only the Yukawa coupling matrix h determines the flavor effect in RH 
neutrino interactions relevant to resonant leptogenesis. 

The RH neutrino mixing effect at resonance can be captured by the one-loop resummed 
Yukawa couplings [18| . Explicitly, in the heavy-neutrino mass eigenbasis, we obtain for the 
resummed Yukawa couplings h^Q, and their CP-conjugates the following values: 


h 




(-2.52-kl.03f) X 10"'^ 
-(2.82-k5.07f) X 10-5 
-(3.63-k4.69f) X 10-5 

-(1.59-k 1.26f) X 10-^ 
(2.71 -k 4.90f) X 10-5 
(3.62-k 4.70f) X 10-5 


(2.39-k l.OOi) X 10-“^ 
(4.37- 3.68i) x 10-5 
(3.77-4.89i) x 10-5 

(1.40 - 1.25i) X 10-^ 
(-l.57-k3.67f) X 10-5 
(-3.78-k4.88f) X 10-5 


(-0.03-k 3.46f) X 10-^ \ 
(-0.08-k 1.49f) X 10-^ 
(0.00 -k 1.84f) X 10-^ / 

(1.02 - 7.03f) X 10-5 \ 
(0.02 - 1.39f) X 10-^ 
(0.00-k 1.84f) X 10-^ / 


( 21 ) 


( 22 ) 


The heavy neutrino decay rates for the Yukawa-mediated two-body processes —> Li4>, 

are calculated with these resummed Yukawa couplings: 


r(iv„ ^ Li^) 


lOvr 


h;„h 


1 


r(Y„ ^ L^icj)^) 


IOtt 




C* 

la * 


(23) 


which are crucial parameters for the calculation of lepton asymmetry, as we will see in the next 
section. 
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4. Resonant leptogenesis with TeV-scale Wr 

As eluded to above, for low scale seesaw models, a simple way to generate enough lepton 
asymmetry is to use resonant leptogenesis [HI EE]. It requires at least two RH neutrinos to 
be nearly degenerate, which is guaranteed in our model by the choice of the texture in the RH 


neutrino mass matrix given by Eq. (12). Thus, in addition to guaranteeing suppression of lepton 
flavor violating processes such as —)• 3e, this choice of RH neutrino texture also provides a 
necessary condition for resonant leptogenesis. As far as the Dirac neutrino mass matrix Mu is 


concerned, we choose the pattern given in Eq. (11) so that it fits neutrino oscillation data. Since 


our focus is on the feasibility of low scale leptogenesis in a realistic model for neutrinos, we do 
not stress here on the naturalness of the texture, although a similar exercise can be performed 
with the fit presented in Section 2.1, The details of leptogenesis in this class of models was 
discussed in [22| and we do not repeat the basic equations here. However, we have extended the 
work of |22] and also present here some revised numerical results due to a computational error 
in the code used there. 

The basic steps in our analysis are as follows: first we calculate the flavored CP asymmetry 


1 


^la — 


riv. 




(24) 


where the partial decay widths are given by Eq. (23) and the total decay width is given by 


= 


I 


[r(iv„ ^ Lict>) + r(iv^ ^ Lf(/)=)] + 2 r(Ar„ ^ iRqm'R) 


(25) 


Note that in LRSM, there is an additional contribution to the total decay width due to the 
three-body decay mediated by Wr: 


T{Na —> IrQrQr) = T{Na —> iRQRq'R) = 






ds 


m 


Nr 


-3ml 


+ 2s-^ 


(S ^ ^Wr^Wr 


(26) 


where Tw^ ~ (g'|j/47r)Mvi/^ is the total decay width of Wr, assuming that all three heavy 


neutrinos are lighter than Wr. The CP asymmetry (24) in our case is determined by the RH 
neutrino texture, specifically, the magnitudes of Mi ^2 and 6M in Eq. (Q as well as the magnitudes 
of the Yukawa couplings, and these parameters are constrained by the neutrino fit. 

The second step in our calculation is to calculate the thermodynamic evolution of the 
normalized heavy neutrino and lepton doublet number densities and respectively, where 
rj^ = IrP and rP = 2m^^('(3)/(7r^z^) is the photon number density, C(x) being the Riemann 
zeta function. For simplicity and for a fair comparison with the results of [2T], here we use the 
flavor-diagonal Boltzmann equations |25j^l 


dr/, 


N 


dz 




-(37-i|(o= + s=), 

^ '/eq 


drjf^^ 

dz 




N 


- 1 


Da - ^Vi 


AL 


Wi 


(27) 

(28) 


where 2 ; = is a dimensionless variable (T being the temperature of the Universe) and 

T/^ = n^q^/rp = z"^K 2 {z)/2Q{3) is the heavy neutrino equilibrium number density, Kn{x) being 

^ Including flavor off-diagonal effects could lead to additional enhancement of the final lepton asymmetry |26[ 127] . 
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Figure 2: Evolution of 'y/{n^^H) [left panel] and 'yK'nl^H) [right panel]. For the latter, we have 
also included the scattering rate of W^l~^ —)• with m^++ = 10 TeV and mwji = 13.1 

TeV. 


the n-th order modified Bessel function of the second kind. The various decay {D^, Da), 
scattering (5^) and washout {Wi) rates appearing in Eqs. (27) and (28) are given by 


Da = 


Da = 


Sa = 


Wl = 


z 

n'yHj\f 

z 

n<Hi^ 

z 

1 




E 




+7^ 


.Sr\ 

OL /’ 


(29) 

(30) 

(31) 




'^lka+ lla + 7i«^ 
k / 






z^K^{z)Kf{z) , 


2C(3) 


(32) 


where = H[z = 1) cz. 17m^^/Mpi is the Hubble parameter at z = 1, assuming only SM 
degrees of freedom in the thermal bath, Mpi = 1.2 x 10^® GeV is the Planck mass and Bia is 
the branching fraction of the RH neutrino decays relevant for the generation of CP asymmetry: 

1 


Bla — 




-[T{Na^Licl})+T{Na^L<i<t>^)] . 


(33) 


The various y’s appearing in Eqs. (29)-(32) represent the reaction rates that involve the RH 
neutrino decays and inverse decays as well as other 2 o 2 scattering processes in the model 


see 


Appendix A). Only the two-body decays of the RH neutrinos involving complex Yukawa 


couplings are responsible for building up the asymmetry, whereas all other processes lead to 
washout effects. 

We present in Figure the evolution of the various collision rates as given in Appendix A 


The left panel of the Figure is relevant to the first Boltzmann equation (27), while the second 


one is relevant to (28). The vertical dashed line shows the critical value z = Zc beyond which the 


sphaleron processes become ineffective. The horizontal dashed line is shown for easy comparison 
with the Hubble rate. Here we have used mwn = 13.1 TeV which is the lowest value of mwn 
we found in this model with an explicit neutrino fit that satisfies the leptogenesis constraints. 

Recently, it was pointed out [3^ that the scattering process W^l~^ —)• W^l'^ mediated by the 
doubly-charged component of the RH triplet might spoil leptogenesis, if and m^++ 
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Figure 3: Various reaction rates '){zc)/{n^'^{zc)H{zc)] as a function of mwji for niN = 585 GeV. 
The green lines in the right panel correspond to the process —)• for various 

masses. 


are relatively small. We have explicitly calculated the washout effect due to this process for 
various values of mwn and m^++. As shown in Figure and also in Figure]^ (right panel), 
this effect is negligible compared with the decay rate and other scattering rates when my/^ and 
m^++ are in the regime of interest. 

In the strong washout regime, the flavored lepton asymmetry is approximately given 
by [Mina 




lOL 


2zKf{zf^'^D^{z) + So,{z) 


Daiz) 


(34) 


where iF®® 


is the effective washout parameter as defined in Eq. (32) and D/{D + S) is the 


effective dilution factor. For the numerical fit given in Section we obtain 


= 

^la 


5.16 X 10“^ 
6.44 X 10“^i 
0 


5.13 X 10“^ 
4.89 X 10“^° 
0 


1.32 X 10“^ 
6.48 X 10“^° 
0 


(35) 


and the total lepton asymmetry is found to be — 2-47 x 10 This is well 

consistent with the observed value of the baryon asymmetry = (6.105i’iQgg®) x 10 “^*^ [29], 
after taking into account the sphaleron conversion rate and the entropy dilution factors. 


4.1. Lower bound of 

Here we present our procedure to derive a lower bound of rnwu which is compatible with 
successful leptogenesis in a LRSM with all RH neutrinos being quasi-degenerate in mass. The 
total lepton asymmetry is given by Eq. (34), which can be rewritten using Eq. (32) as 



^ 2 / \ 1 Da{z) 

4C(3)^ ^^"^^Wi{zf^D^iz) + Sa{z) • 


(36) 


Eor the purpose of deriving a lower bound on mwn, we assume that the dominant washout 
effect comes from the IVR-mediated scattering processes, i.e., Sa ~ Ylk^Ro‘ ~ ^Sjia where 
Sjia = This follows from the fact that gauge interactions are flavor-blind, and 

hence, 7 ^^ has no dependence on lepton flavors. Similarly, when all the RH neutrino masses are 
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quasi-degenerate, we can write Sr = SRa ~ 35rq,. We also assume that D = Ya ~ 
for simplicity. Note that Sr is not summed over lepton flavors while D is. We further define 
and as the CP asymmetry and branching ratio, respectively, without the 3-body decay width 
included in the denominator of Eqs. (24) and (33). For simplicity, we further assume that the 


branching ratios of the decay process Na —>■ Li(p are the same for all the lepton flavors. With all 
the above-mentioned reasonable assumptions, we can approximate the total lepton asymmetry 


in Eq. (36) at the critical temperature as 




4C(3) 

9 z^Ki{zc) 


'^i(^c) ^ ^ [R n 


4 “ SRct^Zc)^ 


4C(3) Sr{zc) {3 + rsrd){3 + rs) 


-e. 


Y 

tot 




Ba i^c) 


Da{Zc) + 3SRa{Zc) 


(37) 

(38) 


where the ratios = D{zc)/S r{zc) and = D{zc)/D{zc) parametrize the relative washout 
strengths of the 2-body decay, WR-mediated decay and VbR-mediated scattering processes. The 
values of and depend on m^r, mw^, and the Yukawa coupling h. Assuming a specific 
value of h and calculating the 2-body decay width simply as T{Na Licj)) = T{Na —>■ L^cj)^) = 

/i^mAr/lGvn we can evaluate these parameters as functions of rriAr and 

.al, 


shows the contour plots of \r]^^(z^)\ = 2.47 x 10 ^ for two different Yukawa 

1-3.5 


Figure 

couplings h = 10“^'® and h = 10“^'®. The red curves correspond to eJot = 1 which is the 
total CP asymmetry in the example fit given in Section]^ Any region outside the red curve is 
incompatible with successful leptogenesis under the assumptions we introduced to obtain it, i.e. 
mjsfi ^ mRf 2 ^ iyins, BJ^ ss BY ~ B^ ~ 1/3, e^t = 1 Ini' specific values of h. For 

h = 10“^'^, we find that the lowest value of mwji allowed for = 1 is 13 TeV at around rriN = 
580 TeV. Note that the example fit we have presented in the previous section has = 13.1 
TeV, rriN = 585 GeV, and the resummed Yukawa coulings of order around 10 ^■®. Thus, 
our example fit, shown as a green dot in Figure]^ is very close to the minimum value of mwH 
obtained here in a simplified manner, which justifies the validity of this approach. In other words, 
we have derived the lower bound in two different ways: (i) in Section]^ by carefully scanning the 
parameter space in order to find an explicit fit with the lowest possible mw^, compatible with 
leptogenesis, and (ii) using the approximate expression in Eq. (|38|). The remarkable agreement 


between these two approaches shows that the simplified expression in Eq. (38) is very effective 


in finding the region in parameter space compatible with leptogenesis and predicting the lower 
bound of rriwn as well as the position in the parameter space where it exists. If we use the same 
expression and take the maximal CP asymmetry allowed in principle, i.e. = 
then we obtain the blue curves Figure]^ With this maximal CP asymmetry, we have found that 
the Yukawa coupling h = 10“^'® gives the lower bound of = 9.9 TeV at rriN = 630 GeV. 

The Yukawa couplings cannot be increased arbitrarily without spoiling the lepton asymmetry, 
since not only the source term due to the two-body decay of the RH neutrinos, but also the 
washout effects due to inverse decay and AL = 2 scattering increase with the Yukawa couplings. 
Similarly, for very small values of the Yukawa couplings, the branching fraction of the two- 
body decay mode becomes comparable or smaller than the three-body decay mode due to Wr 
interactions, and therefore, the washout effect again increases. Thus, successful leptogenesis 
works only in a range of the Yukawa coupling parameter space. This is shown in the three- 
dimensional plot given in Figure where we see that leptogenesis constraints in our model 
require the Yukawa coupling to be 10“®'® < h < 10“^'^ for mwji < 30 TeV. The robustness of 
the lower bound on obtained in Figure]^ can also be verified from Figure]^ 
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Logio[wAr/TeV] 


Figure 4: Contour plots of \r]^^{zc)\ = 2.47 x 10“^ for h = 10“^'^ (dashed lines) and h = 10“^'^ 
(solid lines) with e^ot = 1 (red lines) and = 3 (blue lines). The green dot corresponds to the 
example fit value presented in Section 1^ 


5. Summary 

In this proceedings, we address two issues related to seesaw models for neutrino masses. The 
first one deals with whether the TeV scale can be naturally in the TeV range without fine-tuning 
of parameters. We present a natural TeV scale left-right model which achieves this goal and 
therefore provides a counterexample to the common lore that either the seesaw scale must be 
superheavy or the active-sterile neutrino mixing must be tiny in a UV-complete seesaw model. 
The second issue addresses the important question: Whether in such low scale models, one can 
have successful leptogenesis, and if so, what constraints are implied by this on the mass of the 
RH gauge boson. In an explicit TeV-scale LR model, with an explicit fermion mass fit, we find 
the lower bound to be 13.1 TeV and for generic models in this class of L-R seesaw with enhanced 
neutrino Yukawa couplings compared to the canonical seesaw case and with maximal possible 
CP asymmetry for each flavor, this bounds becomes 9.9 TeV. 

Acknowledgement 

The work of P. S. B. D. is supported by the Lancaster-Manchester-Shefheld Consortium for 
Fundamental Physics under STFC grant ST/L000520/1. The work of C. H. L. and R. N. M. is 
supported by the NSF grant No. PHY-1315155. 


13 





1.0 -3 


(TeV) 


Figure 5: The allowed region of parameter space (shaded region) yielding successful leptogenesis 
in our L-R seesaw model. The vertical gray surfaces show the bound on Yukawa couplings, while 
the horizontal pink surface shows the lower bound on mwn- 


Appendix A. Collision Rates 

The various collision rates for decay and scattering processes in Eqs. (27) and (28) are given 
below: 






Ya = 

it = 
ta = 




NqlLi I NctU^ I NolQ I NocLi I 

iQu^ + iL^Q^ + lL,t + + lL,4> 

N N 

NcLi , Ncu<^ , ^ Nc,Li , 1 

^nIQu^^ ^1LiQ‘= ^ILiu ^1Li 4> ^ll- 

'/eq Veq 


Sr _ NcIr _j_ NoUr NccLr 
'urcIr UrcIr Irur 

la ^.NglR NaUR NadR 
^URdR ' hRdR Irur ’ 

'/eq 

'Siipi I LiL,, 


ha = 

it = 


(AL=2) 

Ilk 

(AL=0) 

Ilk 


'-^kVk 


iLk^k + lLk<pl ^ 


LiL% 

i 


(AT) 

(A.2) 

(A.3) 

(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 


The scattering terms involving two heavy neutrinos in the initial state, e.g. induced by a 
t-channel Wji or eji, and by an s-channel Zji, are not included here since their rates are 
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doubly Boltzmann-suppressed and numerically much smaller than the scattering rates given 
above 1211 ED- The decay rates are explicitly given by 


1li4, 


Nc 


m 


No. 


tt'^z 


Kliz) LicP) + T{Na 


L 




m 


No. 


'k‘^Z 


Ki{z) [r(iVa —)• iRqR^n) + r(iVa —>■ iRQRf^R)^ ■ 


The various collision terms for the 2 •(-)> 2 scattering processes XY 0 AB are defined as 


o.,NY _ 
lAB — 


m 


Ni 




roo 

/ dXy/xKi{zy/x)a^Q{x), 

'^^thr 


(A.9) 

(A.IO) 


(ATI) 


where x = s/m\j^ with the kinematic threshold value Xthr = max[(mx+?T^y)^, {mA+rnB)'^]/ 
and a^^{x) are the relevant reduced cross sections, whose explicit expressions can be found 
in [22]. The reduced cross section of the additional process -A W^l~^ mediated by 

is given by 


a 


WrIr/\ _ 
WrIrY) - 


9Rm%^ 

Stts 


[s — m 


log 


Wr 


+m^++s 


+ 


m ++S 


m 




(s-^Wr. 


2 


- 1 


(A.12) 
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